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ABSTRACT 



The decay B Xs7 is studied at the Next to Leading Order in QCD in a class 
of models containing at least two Higgs doublets and with only one charged Higgs 
boson non-decoupled at low-energy. The two-loop matching condition is calcu- 
lated and it is found to agree with existing results. The complete dependence of 
the Wilson coefficients on the matching scale is given. The size of the Next to 
Leading Order corrections is extensively discussed. Results for branching ratios, 
possible CP asymmetries and lower bounds on the charged Higgs mass are pre- 
sented when the convergence of the perturbative series appears fast enough to 
yield reliable predictions. Regions in the parameter space of these models where 
the Next to Leading Order calculation is still not a good approximation of the 
final result for these observables are singled out. 
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1. Introduction 



It is well known that charged Higgs contributions to the branching ratio for the decay 
B — >• , BR(_B — >• Xs'~i) , decouple very slowly from the Standard Model (SM) one. 
Hence, in the absence of any experimental evidence for a charged Higgs boson, this decay 
may provide a powerful tool to limit the range of unknown parameters in models where 
such a particle is present and other -non standard- contributions to BR(i? — > Xs7) are 
subleading |l|J^. 

Supersymmetric models constitute, perhaps, the best motivated extension of the SM 
where a second Higgs doublet, and therefore a charged Higgs boson if^, is necessary 
for internal consistency [Q. Simpler extensions where only Higgs doublets are added to 
the SM are, however, important on their own right. They are also excellent pedagogical 
tools to understand the subtleties that the Next to Leading Order (NLO) calculation of 
BR(i? — > Xs7) entails, and which may be hidden in the SM results. 

The simplest class of such extensions, with two Higgs doublets, is usually denoted as 
2HDMs. This class contains the well-known Type I and Type II models in which the same 
or the two different Higgs fields couple to up- and down-type quarks. For what concerns 
us here, Multi-Higgs Doublet models can be included in this class, provided only one 
charged Higgs boson remains light enough to be relevant for the process B —>■ Xsj . This 
generalization allows a simultaneous study oi B ^ Xgj in different models, including 
Type I and Type II, by a continuous variation of the (generally complex) charged Higgs 
couplings to fermions. (No tree-level flavour violating neutral couplings are assumed in 
the present paper.) It allows also a more complete investigation of the question whether 
the measurement of BR(i? Xg'-f) closes the possibility of a relatively light not 
embedded in a supersymmetric model 0]. 

At present, a measurement of this decay rate by the CLEO Collaboration is avail- 
able §: _ 

BR(B X,7) = (2.32 ± 0.57 ± 0.35) x 10"^ . (1) 

There exists also a preliminary result by the ALEPH Collaboration with a larger central 
value p: 

BR(B -> X,7) = (3.29 ± 0.71 ± 0.68) x 10"^ . (2) 

Adding statistical and systematical errors in quadrature, one obtains 90% C.L. lower 
and upper limits on BR(5 ^ X,7) of 1.22 x 10"^ and 3.42 x 10"^ from the CLEO 
measurement and 1.68 x 10~^, 4.90 x 10~^ from the ALEPH result. The band of allowed 
values, corresponding to a more conservative estimate of the systematic error, is reported 
by CLEO to be (1.0 - 4.2) x 10"^ §. 

The theoretical situation within the SM, is at the moment far better settled than 
the experimental one. After the original observation that QCD corrections to the decay 
B —>■ Xsj are a substantial fraction of its rate 0, a collective theoretical effort of almost 
a decade has led to the determination of BR(i? —>■ Xs'~f) at the NLO in QCD, and to 



a considerable reduction of the theoretical error. For the Leading Order (LO) calcula- 
tion, several groups contributed to the evaluation of the elements of the the anomalous 
dimension matrix p],^ and provided phenomenological analyses with partial inclusion of 
some NLO contributions (e.g. bremsstrahlung corrections) [pil|, pT|JT^JT3[| . The two-loop 
matching condition, needed for a complete NLO calculation, was first obtained in and 
later confirmed in |T5|JT6|Jr^ , using different techniques. The two-loop corrections to the 
matrix elements were calculated ||T8| and the determination of the 0{a1) elements of the 



anomalous dimension matrix, started already since some time [[191 , ^01 , has been completed 
only recently [^J2^ . In addition, non-perturbative contributions to BR(i? — >■ Xsj) , scal- 
ing like 1/ml |p3|j2^p5| and I /ml p6|, were also computed. The issue of the dependence 



of the branching ratio on scales, first raised at the LO level in [^JT^, was addressed 



for the NLO calculation, and discussed in detail, in [^]. In the SM, uncertainties due 
to sensitivity of BR(_B — > Xsj) on such scales are, undisputedly, rather small and the 
theoretical estimate for this observable suffers mainly from the large experimental errors 
of the input parameters. 

It would be desirable to have similarly precise calculations also in extensions of the 
SM. Thus, more accurate experimental measurements, when available, could provide strin- 
gent constraints on the free parameters of these models. We present here a detailed study 
of B Xg'y at the NLO in QCD, in the class of models specified above, aiming, in partic- 
ular, to an assessment of the reliability of the theoretical calculation. Our results are, in 
general, less optimistic than one could have foreseen. Indeed, we find unexpectedly large 
NLO corrections and scale dependences in the Higgs contributions to BR(i? —>■ Xsj) , ir- 
respectively of the value of the charged Higgs couplings to fermions. This feature remains 
undetected in Type II models, where the SM contribution to BR(i? Xg'-f) is always 
larger than, and in phase with, the Higgs contributions. It can, however, produce un- 
acceptably large scale uncertainties for certain ranges of these couplings and, at times, 
completely ill-defined results. We single out combinations of couplings and of values of 
the charged Higgs mass, m^^, for which the branching ratio can be reliably predicted. In 
this case, a comparison between theoretical and experimental results for BR(i? — ^ Xg'y) 
allows to conclude that values of uIh = 0{M^-) can be excluded, in general, only in 
Type II models, but are otherwise allowed. 



Previous LO analyses had dealt with Type I and Type II models P, p5| , pU| , pU | and 



the generalized class of models considered here pT|j3^ . LO calculations are known to 
have large scale uncertainties. Consequently, they are not particularly good arenas to 
distinguish the quantitative differences between QCD corrections to charged Higgs and SM 
contributions. These differences stand out clearly at the NLO level. Two NLO calculations 



have been performed recently |[T7|J33[| . They deal with Type I and Type II models and 
the issue of scale uncertainty is addressed in for Type II models. The fact that the 
NLO corrections to the Higgs contributions are large is, therefore, understandably missed, 
since the branching ratio is dominated by the SM contribution. From the technical point 
of view, our calculation of the two-loop matching condition agrees with that reported 



in p!7| , p3| . We disagree, however, with the analytic dependence of BR(i? -^^7) on the 
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the matching scale given in the hterature for the SM 



PI 



and for 2HDMs 



|T7| 



The reminder of this paper is organized as follows. In Sect. ^ we define the class 
of 2HDMs studied. In Sect. ^ we outline the main steps of the calculation; to keep this 
section readable, we relegate some parts of our results to appendices. In Sect, ^we present 
our phenomenological analysis of the decay B —>■ Xg'j in the class of models considered. 
After giving our branching ratio prediction for the SM, we discuss, in [4.1| , the size of the 



NLO corrections to the charged Higgs contributions at the amplitude level. In Sec. |4.2| , 
we give results for BR(i? Xg'y) , and discuss their reliability, for various ranges of real 
couplings of the charged Higgs to fermions. NLO branching ratios for Type I and Type II 
models are in particular studied in Sec. O, where lower bounds for m^, are given within 
Type II models. In Sec. we consider the decay B — > when complex couplings 
are involved and we investigate the impact of their phases on CP rate asymmetries of 
B —>■ Xsj . We give also specific examples of couplings for which the theoretical prediction 
of the branching ratio is reliable and compatible with the CLEO measurement, even for 
rriH = 0{Mnr). Finally, our concluding remarks are in Sect. ^ 



2. Two Higgs Doublet Models 

Models with more than one Higgs doublet have generically a Yukawa Lagrangian of the 
form: 

- £ = /ij 01 d'n, + 02 u'j,^ + 03 e'«^. + h.c. , (3) 

where g^, 0i, {i = 1,2,3) are SU(2) doublets (0i = ia^cj)*); u'^, rf'^, and e'^ are SU(2) 
singlets and h'^, /i", and U denote 3x3 Yukawa matrices. Gauge invariance imposes 
the value of hypercharge Y{(j)i) = 1/2. Suitable discrete symmetries are usually invoked 
to forbid additional terms like hfj q'^^ 02 c^'ijj which would induce flavour changing neutral 
couplings at the tree-level. Indeed, to avoid them altogether, it is sufficient to impose 



that no more than one Higgs doublet couples to the same right-handed field p4 | 



When only two Higgs doublets are present, 0i and 02, it is in general 03 = 0i (or 
03 = 02)- Nevertheless, for the sake of generality, we leave the symbol 03 distinct from the 
other two. We indeed include in our discussion, as "effective" 2HDMs, also models where 
an n-number of sequential Higgs doublets is present and we assume that the additional 
charged Higg bosons other than the lightest one become heavy enough to decouple from 
our problem. 

After rotating the fermion fields from the current eigenstate to the mass eigenstate 
basis, the charged Higgs component in eq. @) becomes: 

~ ~7~S~Z "^LjVji 01+ - "^0' UR.Vij dij 02^ + "^0' T^Li eRi 03+ + h.c. , (4) 

< (pl > < ^ ^ ^3 ^ 

where Vji are elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix. A further 
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rotation of the charged Higgs fields to their physical basis through a unitary matrix f/, 



U 



(5) 



yields the following Yukawa interaction for H^: 



XuLjVji dm + 



YuRiVijdLi + 



mii 



ZVLieRAH+ + h.c.. (6) 



All fields Tij^ are supposed to be heavy enough to become irrelevant for phenomenology 
at the electroweak scale. In (j^), the symbols X and Y are defined in terms of elements of 
the matrix U 



35.36.37.31.3 



X 



u 



12 



Y 



U22 
U21 ■ 

-U32/US1 if 



(7) 



The symbol Z has a similar definition, i.e. Z = if 03 7^ 0i, 4>2 or coincides with 

X (— Y) if 03 = (pi (03 = 02)- Notice that X, Y, and Z are in general complex numbers 
and therefore potential sources of CP violating effects I^S],^. Their values are only very 
loosely constrained by the requirement of perturbativity and low-energy processes such 
as the B-B mixing ||32| . 



(8) 



When only two doublets are considered, the diagonalization matrix f/ is a 2 x 2 
orthogonal matrix 

/ cos /3 — sin /3 \ 
y sin f3 cos l3 J 

Although both doublets are present in the theory, one still has the freedom of selecting 
01 = 02 in eq. (^). This choice gives rise to the 2HDM of Type I, to be distinguished from 
the Type II in which both doublets contribute to the Yukawa interactions |^ . It is easy 
to see that in these two cases the couplings X and Y are real and given by: 



X = -cot/3 , 
X = tan/5 - 



Y 
Y 



cot/5 
cot/5 



(Type I) 
(Type II) 



(9) 



Note that the coupling of the Goldstone boson 0+ to matter fields is independent of 
the number n of Higgs doublets considered, and always equal to: 

The calculation of the B — > Xsj decay rate in these models, therefore, is modified simply 
by the addition of charged Higgs contributions to the usual SM one. The results can 
be described in terms of m^^ and the two complex parameters X and Y. The presence 
of phases in the couplings X and Y allows for CP asymmetries in the decay rate for 
B^Xs^. 
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3. The Calculation 



We use the framework of an effective low-energy theory with five quarks, obtained by 
integrating out the heavy degrees of freedom, which in the present case are the t-quark, 
the VT-boson and the charged Higgs boson. As in the SM calculations, we only take into 
account operators up to dimension six and we put rris = 0. In this approximation the 
effective Hamiltonian relevant for radiative i?-decays (with |Ai?| = \AS\ = 1) 

^eff = -^V,:v,,J2Qift)0,i^) (11) 
v2 

consists precisely of the same operators Oi(/i) used in the SM case, weighted by the Wilson 
coefficients Cj(/i). They read: 



(12) 



05= {sL7,iuipbL)E,{qrrYq), 06= {sLl,l.%T%L)E,{qrYYTaq) 
07 = mtifi) (sLa^^'-bR) F^, , 08 = if^ m,(/i) {-SLa>^-T%R) G% , 

where (a = 1, 8) are SU{3) colour generators; Qs and e, the strong and electromagnetic 
coupling constants. In eq. (0), mb{fj,) is the running 6-quark mass in the MS scheme 
at the renormalization scale /i. Henceforth, mq{^) and iriq denote MS running and pole 
masses, respectively. To first order in a^, these masses are related through: 

m» = ™il + ^ln4-i^^*)). (13) 

y TT 6 -K ) 

Note that the equations of motion have been used when writing down the list of operators 
in eq. (|12|). This is sufficient since we are interested in on-shell matrix elements |^T| and 
since we choose to perform the matching by comparing on-shell amplitudes obtained in 
the effective theory with the corresponding amplitudes in the full theory. The reader 
who is interested in doing off-shell matching (and therefore working in a (larger) off-shell 
operator basis) is referred to ref. |]T7 . 



Working to NLO precision means that one is resumming all the terms of the form 
a^inib) \v^{mb/ M), as well as asimb) {a^imb) ln"(mf,/M)), where M stands for one of 
the heavy masses M^, rrit ot rriH- This is achieved by performing the following 3 steps: 

1) One matches the full standard model theory with the effective theory at a scale ^yy, of 
order M. At this scale, the matrix elements of the operators in the effective theory lead 
to the same logarithms as the calculation in the full theory. Consequently, the Wilson 
coefficients Ci{fiw) only pick up "small" QCD corrections, which can be calculated in 
fixed order perturbation theory. In the NLO program, the matching has to be worked out 
at the 0(aq) level. 
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2) The evolution of these Wilson coefficients from /i = /x^ down to fi = fib, where fib is of 
the order of m^, is obtained by solving the appropriate Renormalization Group Equation 
(RGE). As the matrix elements of the operators evaluated at the low scale /ib are free of 
large logarithms, the latter are contained in resummed form in the Wilson coefficients. 
For a NLO calculation, this step has to be performed using the anomalous dimension 



matrix up to order 



3) The corrections to the matrix elements of the operators (s7|Ci(/i)|6) at the scale = fib 
have to be calculated to order as precision. 

The charged Higgs boson enters the NLO calculation only via step 1). The Higgs boson 



contribution to the matching condition is obtained in the same way as the SM one |T5 
therefore, we do not repeat any technical details. A general remark, however, is in order. 
In the procedure described above all heavy particles (W^-boson, t-quark and charged 
Higgs) are integrated out simultaneously at the scale fiw In the context of a NLO 
calculation this should be a reasonable approximation provided is of the same order 
of magnitude as or rrit. 

Before giving the results of the three steps listed above, we should briefly mention 
that instead of the original Wilson coefficients Ci{fi) is it convenient to use certain linear 
combinations of them, the so-called "effective Wilson" coefficients C^^{fi) introduced 
in 



13,21 



C;^{fi) = Qifi), (^ = 1,...,6), 

Cfifi) = CM + Y.y,Q{fi), Cf{fi) = CM + Y.z^C^{^^), (14) 



i=l 



where Ui and Zi are defined in such a way that the leading oder matrix elements {s'^\Oi\h) 
and {sg\Oi\b) {i = 1,6) are absorbed in the leading orders terms in CY^{fi) and CQ^{fi). 
The explicit values of {y,} and {zi}^ = (0, 0, -|, -|, -f , -f ), ^ = (0, 0, 1, -|, 20, -f ) 
were obtained in ref. |2T[] in the MS scheme with fully anticommuting 75 (also used in 
the present paper). 



3.1. NLO Wilson coefficients at the matching scale fi^: {fiw) 

To give the results for the effective Wilson coefficients C^^ at the matching scale Hw in a 
compact form, we write: 

cr{^^^) = c^^-'M + ^^cl^-'M . (15) 

The LO Wilson coefficients at this scale are well known |^2|,^. We decompose them in 
such a way to render explicit their dependence on the couplings X and Y: 
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C°'^'^(M = (z = 1,3,4,5,6) 

Cs'^^ifJ'w) = C's.sM + lYpCgyy + (XY*) Cg . (16) 

The coefficients Cj^^j^j{nw) and Clgj^ji^Hw) are functions of x = m'^/M'^, while Cyj(/ii^,-) 
and C^jifiw) (j = YY, XY) are functions of ?/ = mf/m'jj; their exphcit forms are given in 
Appendix |^. Note that there is no explicit dependence of the matching scale fiw in these 
functions. Whether there is an implicit /Xvi^-dependence via the t-quark mass depends on 
the precise definition of this mass which has to be specified when going beyond leading 
logarithms. If one chooses for example to work with ml{fi^), then there is such an implicit 
^vi'^dependence of the lowest order Wilson coefficient; in contrast, when working with the 
pole mass rrit there is no such dependence. We choose to express our (NLO) results in 
terms of the pole mass m^. 

The NLO pieces Cl''^^{^„) of the Wilson coefficients have an explicit dependence on 
the matching scale fi^ and for i = 7,8 they also explicitly depend on the actual definition 
of the t-quark mass. Initially, when the heavy particles are integrated out, it is convenient 
to work out the matching conditions Cl'^^^^fiw) for i = 7, 8 in terms of rrit{^yv)- Using 
eq. ([T3|), it is then straightforward to get the corresponding result expressed in terms of 
the pole mass irit. As in the LO case we give them in a form where the dependence of the 
couplings X and Y is explicit: 

Yl,eff/ N 1 r- , 1 l^W 



Cr^'M = 15 + 6 In 



M2 

2 



Ci'^"(M = i?o + §ln^ + |Y|^EH 

Cl'^^'M = (2 = 2,3,5,6) 

Cg^^it^w) = Cl'J^j{flw) + |Yp C8,yy(/^w) + (X Y*) Cg'^^ (/ivi/) , (17) 

where for i = 7,8 the three terms on the right hand side can be written in the form 

Cl'stM = iy,sM + M,,,,ln^ + T,,M(ln^-^) 

CI:^^yM = 1^.,.. + M,,,,ln^ fln^-^) 

CI'^yM = W^.,x. + M,,, In % + T,,, fin ^ - ^] . (18) 

Note that in eq. (|18D the Wij- and the Mij (j = SM,XY,YY) terms would be the full 
result when working in terms of mt{fiw)- The Tij terms result when expressing fritifJ'w) 
in terms of the pole mass rrit in the corresponding lowest order coefficients. Thus, for 
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i = 7, 8, the Tj ,,— quantities are 



n5M = 8x , T,, = 8y — , (j=XY,YY). (19) 



Notice that if one worked with the running t-quark mass mtifJ't) normahzed at the scale 
Ht instead of the pole mass rrit, the third terms on the right hand sides of eqs. (^) would 
have to be replaced by Tj j ln(/ij//i^) (j = SM, YY, XY). The functions Wtj, Mij and Tj j 
(j = SM,XY, YY), together with Eq and E^, are listed in Appendix ^. Our results for 
CY''^^{fiw) and Cg'^^{fiw) agree with those in refs. and [0, when taking into account 
that the latter results are expressed in terms of mt{[iw)- We disagree, however, with 
the form of the coefficients C\'^^{^,y^) and C\'^^ {^.w)-, whose matching scale dependence 



is forgotten in ||27| , p!7[| . As it will be seen later on, this neglect has consequences in the 



determination of the matching scale which yields the lowest estimate of BR(i? Xs'~f) . 



3.2. NLO Wilson coefficients at the low— scale /i^: Cl^{iib) 

The evolution from the matching scale down to the low-energy scale /i^ is described 
by the RGE: 

l^j-^Cf{^,) = Cf{^,)^;f{^,). (20) 

The initial conditions Cl'^{nw) for this equation are given in Section pTT| , and the anoma- 
lous dimension matrix ■y^^ is given in Appendix up to order a|, which is the precision 



needed for a NLO calculation. The solution of eq. (|2(]|) , obtained through the procedure 
described in [|^, yields for the coefficient 



cf{^,) = c^-'{^^,) + ^cl-'{^^,), (21) 



the LO term: 



1=1 

and the NLO one: 

r^l.cff/ \ 39 1 fj X 8 / 37 39N 1 fj 

/297664 16 7164416 n 256868 37 6698884 39\ n , 

+ 7723 7723 -I 7723 7723 c; ' (u ) 

V 14283 ' 357075 ' 14283 ' 357075 ' / ^ ' 

37208 / 39 i6\ ^0 eflf/ \ 

+ i76rv''""''")^°^^^-) 

+ E (e^r]Cl'^^(/iH') + (/. + M) C°2'"''M + hvCl'^^'ifi^)) v''' , (23) 

i=l 
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The symbol r] is defined as 77 = as{fiw)/c(s{f''b)] the vectors Oj, hi, Cj, /«, fcj and k are 
listed in Appendix |y. Notice that eq. ( P^D can be used in this form for all models in 
which the same set of coefficients Ci{fiw) are non vanishing. It is more general than the 
corresponding eq. (21) given in which can be used only when the matching scale fiy^ 
is fixed at the value Mw The use made of that equation in |2^, where the dependence 



of BR(i? — s> Xs7) on the matching scale is studied, is inappropriate. 

As far as the other Wilson coefficients are concerned, they are only needed to LO 
precision in the complete NLO analysis of BR(i? X^'y) . In this precision, the Wil- 
son coefficients of the four-Fermi operators {i = 1, ...,6) are the same as in the SM. As 
the coefficients C^^{fib), ...,CQ^{iJ,b) are numerically much smaller than C2^{fj,b), we ne- 
glect contributions proportional to these small Wilson coefficients in the amplitude for 
BR{B — >• Xs'j) and list here only the LO expressions for Ci^{fih), C2^{fj,b), and Cg^^fib): 



c^s'^'^M = v^^c's''''M + Y.Kv<c','^''M, (24) 

i=l 



When discussing the contributions due to the charged Higgs boson, it is convenient to 
spht the Wilson coefficients at the scale fib into the contributions C^^j^{fib), Ci'^Yif^b), 
and Cl^yifih)- 

Cf{l^b) = Cf,M + |Ypq^,(^,) + (X Y*) C5,(/.,) . (25) 

As the solution ( p3D of the renormalization group equation ( pO[ ) is linear in the initial 
conditions, this splitting is given in an obvious way in terms of the corresponding splitting 
at the matching scale /Xvk, presented in Section pj 



When calculating NLO results in the numerical analyses in Section ^, we use the 
NLO expression for the strong coupling constant: 

with 

where /?o = y and Pi = However, for LO results we always use the LO expression 
for as(/i), i.e., /?i is put to zero in eq. (|2 



3.3. Branching ratio for BR(5 -^sT) 

We first give the formulas for the quark decay b Xg'j and discuss the meson decay 
B — >■ later. In a NLO calculation, b Xs'y involves the subprocesses b ^ s'y 
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(including virtual corrections) and b s'jg, 
amplitude for the first can be written as 



i.e., the gluon bremsstrahlung process. The 



A{b S7) 



where the reduced amplitude D is 



D = C°'^^(^,) + 



'72 



V;,Vt,D{sj\Or\b) 



tree j 



An 



The symbol V{fib) is defined as: 



{2i 



(29) 



(30) 



In writing eq. (^) we directly converted the running mass factor rnb{fib), which appears 
in the definition of the operator Oj in eq. ([1^ , into the pole mass mi, by making use of 
eq. ([13|) . This conversion is absorbed into the function V{fib) and consequently the symbol 
{s'y\Oi\b)tree IS the tree-level matrix element of the operator Oi, where the running mass 
factor mh{fib) is understood to be replaced by the pole mass nit,. In the previous literature, 
this procedure was done in two step. First, mb{[Jth) was expressed in terms of rrTbimb) and 
then in turn rfibimb) was converted into the pole mass m;,. This last step brought into the 



game the quantity F (see e.g. ref. [0). The elements 7^7'^^ of the anomalous dimension 
matrix and the virtual correction functions in (|30|) are given in Appendix and 0, 
respectively. Note, that some parts of the bremsstrahlung contributions associated with 



07 are effectively transferred to r^ as detailed in ref. |18 . 



A splitting analogous to that in ( |25| ) holds also for the reduced amplitude D: 

D = DsM + \Y?Dyy + {XY*)D^y . (31) 
as well as for the function V{^b). 

From A{b — > S7) in eq. (|28|) the decay width V{b 57) is easily obtained to be 



T{b 57) 



32^ 



\V^\Vtb\^ae^ml\D\^ . 



(32) 



In the numerics we discard those terms in \D\'^ which are explicitly of order a|. Note, 
however, that there are implicit higher order terms which are retained. Such terms arise for 
example because we evaluate the quantity rj = {as{fiw)/o:s{fib)) using the NLO expression 
for as also in those Wilson coefficients which are needed only in LO precision. 

To obtain the inclusive rate for b Xg'j consistently at the NLO level, we have 
to take into account the bremsstrahlung contributions [llO| , ^ . The corresponding decay 
width r(6 ^ S'jg) is of the form 



T{b s-fg) 



32^ 



Op 



(33) 
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where A is 



A 



71 



0,cS 



0,eff 



(/^b)]* fij} 



(34) 



As in the virtual contributions we put C^"'™ = for z = 3, 6. In contrast to ref. |2T|, 
we do not introduce a cut-off when the photon gets soft. In order to cancel the infrared 



singularity, which appears in this case, we include as in ref. |]T8| the virtual photonic 
correction to the process b sg, which we absorb into the quantity /gg (see eq. (B.9) 
in ref. [|l8l). Note also, that in our approach the term fn is already absorbed into the 
function rj. The non-vanishing fij terms, which have to be taken into account explicitly 
in eq. ( ^if ) , are listed in Appendix |^. 



In order to get the decay width for the meson decay B 
the non-perturbative corrections which scale like l/i 
l/ml [^. The decay width T{B Xs'y) then reads: 

^.7) 



Xg'-f we take into account 
and those which scale like 



r B 



3271^ 



Vts^tb\ Olem 



\Dr + A + 



^|C°'^«^(/x,)p + 



^NP 

mi 



0,cS/ 



1 



yO,cS 




(35) 



Since we work in the new operator basis, a contribution to the term in 1/ comes not only 



from O2, as it was incorrectly assumed in ||T^, but also from Oi. The non-perturbative 
2|] and (5f ^ ^ in (||) are: 



quantities 5. 



NP 
7 



^.P 



(36) 



Ai 9A2 ^jvp _ _'^2 

where Ai and A2 parametrize the kinetic energy of the 6-quark and the chromomagnetic 
interactions, respectively. Their values are: Ai = — 0.5 GeV^ and A2 = — 0.12GeV^. Of 
these two parameters, the first has larger uncertainties that the second one. As it will 
appear from the final formula for BR(i? Xg'^) , the overall Ai dependence cancels to a 
large extent. 



The branching ratio is then obtained as 



BR(5^X,7) =^^t^^^BR 



SL y 



(37) 



where BRg^ is the measured semileptonic branching ratio and the semileptonic decay 
width is given by: 



6i 



m^ 

Z3- (3^ 



mt 



The phase space function g{z) and the (approximated) QCD-radiation function f{z) 
in (|38|) are: 



Iz + 8z^ - z^ - I2z^ In z . 



(39) 
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and the non-perturbative correction reads [p^ , p5[] : 



Ai 
2 



3A2 
2 



4- 



(40) 



In eq. ( pq ) is the renormahzation scale relevant to the semileptonic process, which 
is a priori different from the renormahzation scale fib in the radiative decay, as stressed 

m 



in [^]. However, as pointed out in ref. 
more conservative. We therefore use Jlj, 
the evaluation of eq. (P7D, we do not expand l/F 



the identification Jlf^ = fib turns out to be 
fib in the numerical analysis. Furthermore, in 
SL in powers of ag- 



4. Results 



We discuss in this section results and theoretical uncertainties of the NLO calculation of 
BR(i? Xg'y) in 2HDMs. As it will be shown, these uncertainties can be very large, in 
contrast with what is found in the SM. For reference, we give our SM result: 



BR(i?^X,7) = (3.57 ±°:?^ (/ib) ±ro8 M ^olr (param)) x IQ-'. (41) 



The central value 3.57 x 10~^ is obtained for fib = 4.8 GeV, fiy^, = and the central val- 
ues of the input parameters listed in Table |l| of Appendix ^ The low-scale variation in the 
interval [2.4, 9.6] GeV gives the maximum value for BR(i? — ^ ^^7) at /if, = 4.2 GeV and 
the minimum at fib = 9.6 GeV. The matching scale dependence of the branching ratio is 
monotonically decreasing for increasing /lyy-. from the central value in (^Tf) , BR.{B — ^ ^^7) 
is reduced by 2% at fi^ = rrit. The value 3.57 x 10~^ in eq. (p]) reduces to 3.46 x 10~^ 
when the factor l/F^^ in eq. (|37|) is expanded in ag- 



Our results for the branching ratio BR(i? — ^ -^^^7) in the class of 2HDMs considered, 
is parametrized in terms of {X, Y, mj:^}. We limit the range of this three-dimensional 
parameter space as follows. We fix Y to small real values of 0{1) and scan X in the 
complex plane without ever violating the perturbativity requirement on X and Y discussed 
p^ . As explained in the value of itIh in 2HDMs can be as low as the LEP I lower 



m 



bound of 45 GeV. We do not strictly apply cuts on {X, Y, m^} due to the measurement of 
Rb or of other processes like B-B mixing with virtual exchange of charged Higgs bosons, 
which could remove some corners of the parameter space considered. Our interest is more 
to show some theoretical features of the NLO calculation of BR(i? — > Xs7) in regions as 
wide as possible and to discuss constraints due to the decay B Xg'y itself. For the 
numerical evaluations, unless otherwise specified, we use the central value of the input 
parameters listed in Table |l| in Appendix |^. The values of the matching and low-scale, 
chosen respectively in the intervals [M^y, Max(mt, m^^^)] and [2.4, 9.6] GeV, are explicitly 
given for each result presented. 

In order to investigate the reliability of the NLO prediction for the branching ratio, 
we study the term |-Dp, which dominates the expression within curly brackets in (35). 
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Notice that \D\^ encapsulates all the matching scale dependence and the bulk of the 
dependence on the low-scale. At the NLO level, the procedure of squaring D dictated by 
perturbation theory, in which terms of are omitted, gives: 

\D\^ = \C^^'''{^^b)?[l + 2Re[AD)] , (42) 

where AD is defined as (see eq. (P^)): 



If |AZ}| is not small, the formally Next to Next to Leading Order (NNLO) term |ADp, 
dropped in eq. (^2|), is numerically relevant. Its omission can lead to branching ratios 
with a large dependence on yU;,, or in extreme cases, even to negative values for the 
branching ratio. In these situations, the truncation of the perturbative expansion of D 
at the NLO level is certainly not justified. As we shall see, the size of |AD| depends 
crucially on the values of X, Y and m^. We split our analysis as follows. In Sect. we 
disentangle the effect of the couplings X and Y by studying the NLO contributions to the 
Higgs components Dxy and Dyy of the reduced amplitude D. In Sect. we illustrate 
predictions for the branching ratio for real X and Y couplings, giving particular emphasis 
to cases in which the results are highly unstable or altogether unacceptable. In Sect. ^]3| 
we discuss in detail Type II models, for which the NLO corrections are under control, 
and Type I models for which the reliability of the theoretical prediction depends strongly 
on the point of parameter space considered. Finally, in Sect. [4.4| , we study branching 
ratios and CP asymmetries in the presence of complex couplings, and outline regions of 
parameter space where these predictions are trustworthy. 



4.1. Amplitude: X— Y independent analysis 

Before analyzing the size of the various corrections within 2HDMs, it is worth reviewing 



the situation in the SM. There, as it is was first observed in |21], the NLO contribution to 



the matching condition tends to cancel the contribution due to the NLO evolution from 
the matching scale to the low-energy scale /if, (see upper frame in Fig. |l]). Indeed, 
which is 12% of C^'g^ at is practically negligible for G [2.4, 9.6] GeV. The 
dominant NLO effect is therefore due to K,A/(yUfe) which, in the same range of /i^, is at 
most 20% of the leading term Cj'jfiifib) and vanishes exactly at fi^ ~ 3.5 GeV. These 
features are illustrated in Fig. ^ The size of VsMif^b) results from cancellations among 



the individual contributions of the four operators Oi, O2, C7 and Og retained in eq. (|30D . 
The details of the above described cancellations and, more generally, of the distribution 
of the complete NLO correction among the individual terms in eq. (pOj), are specific to 



the MS scheme (with anticommuting 75) and could be altered by a different choice of 
scheme. We remind that only the complete -Dsm term as well as the final BR(i? — ^ Xg'j) 
are scheme-independent, up to higher order. 
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Figure 1: RGE evolution of the Wilson coefficients Cf^, at the LO (with NLO as) (dashed lines) and 
at the NLO (solid lines), for [l^ = Myy. The upper frame shows the SM coefficients, the lower one the 
coefficients Cy^^y, Cy^^y, for = 100 GeV. The needed input parameters are fixed at their central 
values listed in Table |^ in Appendix ^. 
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Figure 2: Low-scale dependence of the term D, (solid line), of the LO coefficient C^''^^ {/J'b) (with 
NLO as) (dashed line), and of Cj^i^Hh) (dotted lines), for = M^. For D xy and -Dyy, the value 
mu = 100 GeV was used. The needed input parameters are fixed at their central values listed in Table 
in Appendix 



One is naturally led to ask whether the cancellations among different sets of NLO 
corrections observed in the SM are spoiled in 2HDMs. To address this issue, we investigate 
the building blocks Dxy and Dyy defined in eq. (|3T|), for a specific value of m„. In Fig. 0, 
we show the real parts of Dxy and Dyy for = 100 GeV, together with the real part 
of DsM- Even when considering real X and Y couplings, imaginary parts to D come from 
the absorptive terms of the loop corrections in (|30|). As it can be seen from eq. (^21), the 



imaginary parts of Dxy, Dyy and Dsm do not contribute to BR(i? — > Xs'y) at this order 
of perturbation theory if X and Y are real. In the remainder of this Section, therefore, 
any reference to these components is understood as a reference to their real parts. 

Unlike in the SM, the NLO corrections coming from Cj^xvif^b) and VxyifJ'b), have 
roughly the same size 
correction 



20%) and the same sign, for = 100 GeV. The combined 



AD, 



D, 



An 



(44) 



amounts to the considerable values of —43% to —36% when varying /i^ G [2.4, 9.6] GeV. 
(We warn here that all components of the LO Wilson coefficient Cj''^^{fib) discussed in 
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Figure 3: Low-scale dependence of the term D, (solid lines), of the LO coefficient Cj''^^ [fif)) (with NLO 
as) (dashed hnes), and of Cf^ (fib) (dotted hnes), for rriu = 500 GeV. The matching scale fiy^ = 
is used in the upper frame and to yU^^ = TTljj in the lower one. The needed input parameters are fixed 
at their central values listed in Table |^ in Appendix ^. 
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this Section and plotted in Figs. |l| and ^ are evaluated with NLO ag.) Differently than in 
the SM, there is no scale fib in the range considered, at which the LO and NLO prediction 
for DxY coincide. Similar results hold for Dyy 

The change in matching scale from fi^^ = (value used for Fig. ^ to, say, = ^h, 
is practically inconsequential for itIh = 100 GeV. It becomes very relevant for large values 
of mil, since it crucially affects the size of the NLO corrections. For = 500 GeV, 
as shown in Fig. ^ the correction ADxy reduces to —37%, —31% when varying /if, in 
the usual interval [2.4,9.6] GeV, whereas values of —51%, —45% are obtained for /i,^ = 
Mn^. Notice that, for = 500 GeV, \Dxy\ 

is roughly ten times smaller than |Z)gjv/| 
and Dyy completely negligible. When using the same matching scale fiw — 500 GeV 
for the SM contribution, the cancellation between the NLO correction to the matching 
condition and the NLO evolution of the Wilson coefficient Cy'^f does not occur anymore. 

Large cancellations are instead observed between VsAiif^b) and Cj'sfj{fib)- The complete 
correction \ ADsm\ ranges still between 19% and 2% for fib G [2.4,9.6] GeV, but the point 
where NLO and LO predictions for Dsm coincide is pushed to the higher end of fib- These 
results imply that for heavy and large enough X and Y couplings (to lift at least Dxy 
to be of 0{Dgnj)), a choice of of 0{mjj), instead than (^(M,^), minimizes the size of 
the NLO corrections. 

We draw attention to the fact that the sensitivity of the reduced amplitude D to 
rriH is weaker than that of the Wilson coefficient C^^^^fiw) at the matching scale. It is 
interesting to see that the coefficients Cj^j^^ and Cj^y are almost identical at the matching 
scale fiw for = 100 GeV. This feature is clearly visible in Fig. |l|, where the matching 
scale fiw = is used. (It remains true for all fiyy G [M^,mt].) We observe in Fig. |1] 
that the RGE flow from fj,^ to fj,b introduces a large gap between Cj^j^{nb) and Cj^y{fib), 
which is then somewhat widened at the level of the reduced amplitude D by the inclusion 
of the l^(/ib)-contribution. At fib G [2.4, 9.6] GeV, Dxy is about a factor of three smaller 
than Dgj^i. 

When going to physical observables such as BR{B — > X^'j) , the individual building 
blocks of D, discussed above in a coupling-independent way, are weighted according to 
the values of the couplings X and Y. 

4.2. Branching Ratio: Real Couplings 

The almost perfect flatness of Dyy and Dxy shown in Figs. ^ and ^, should not lead 
to the conclusion that the NLO prediction for BR(i? — > X^'-f) is well-behaved. It was 
explicitly shown in [^, by expanding the Wilson coefficients around fib = rrib, that the 
dominant scale dependence of the form as{mb)lia{fib/fnb) is cancelled in the complete 
NLO expression for D in eq. ( [29| ) and consequently also in |-Dp. If the NLO corrections 
are large enough to reduce substantially the magnitude of the LO term, |_Dp becomes 
sensitive to higher order dependence on fib of the form aKnib) IvP^fib/nib) {p = 1, 2). The 
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Figure 4: Branching Ratio as a function of the couphng X, for Y = 1 and TTljj = 100 GeV. The upper 
(lower) frame shows the LO (NLO) result for /it = 4.8 GeV (solid lines), /ib = 2.4 GeV (dashed lines) 
and fib = 9.6 GeV (dash-dotted hnes), and matching scale fi^ = "^h = 100 GeV. The needed input 
parameters are fixed at their central values listed in Table |l| in Appendix ^. Superimposed is the range 
of values allowed by the CLEO measurement. 
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flatness of D seems to indicate tliat tlie omitted NNLO term |ADp would cancel to a 
large extent this remaining fib dependence. Nevertheless, whether the omitted |ADp in 
eq. ( ^2]) is the bulk of the NNLO corrections can only be decided when a complete NNLO 
calculation is at hand. 

It is clear that the reliability of the NLO prediction for the branching ratio, which 
is linked to the size of AD, depends on the values of X and Y. For a given rriH, it is 
possible to choose these couplings in such a way that the reduced amplitude D is dom- 
inated, for example, by the Higgs contribution Dxy For such points in the parameter 
space {X, Y,mjj}, the size of the NLO corrections to the branching ratio is then roughly 
2Re{ADxY), i-e. about —80% for ttih = 100 GeV. In this case, the NLO corrections 
substantially reduce the leading order prediction. As expected, the resulting scale depen- 
dence of the branching ratio BR(i? — Xg'y) is large, viz. about 40%. For Y = 1 and the 
same value of m^, Dxy completely dominates the branching ratio for |X| > 20, outside 
the range shown in Figs. ^ The prediction for BK{B —>■ Xg'-f) is then, however, far above 
the band of values allowed by the CLEO measurement. 

Moreover, it is possible to choose X and Y in such a way that Re(AD) < —50%. 
As eq. ( ^21) shows, this choice leads to a negative NLO prediction for BR{B — > Xs'j) . 
Needless to say, in such a situation higher order corrections to the NLO calculation are 
mandatory to obtain sensible results. In our representative case of = 100 GeV and 
Y = 1 illustrated in Fig. ^ the range of X corresponding to a negative branching ratio is 
roughly — 10<X<— 2. A comparison with the upper frame of this Figure shows that, 
given the very large scale dependence of the LO calculation, one could have guessed the 
pathological situation which is encountered at the NLO level. 

Reliable predictions for the branching ratio with a mild scale dependence are obtained 
only for X > — 1 (see lower frame of Fig. For these values, the SM contribution 
dominates, but the Higgs contribution proportional to X Y* is still large enough to produce 
a sharp raise of BR{B Xg'j) when X increases from —1 to 2. Notice that for Y = 1 
the values X = — 1 and X = 1 correspond respectively to the ordinary 2IIDM of Type 
I and Type II with tan/5 = 1. Type I models, however, are not always so stable and 
well-behaved as in the case described above. For Y = 1 and X = — 1 (corresponding to 
tan/? = 1), it is enough to lower to 45 GeV to find a low-scale dependence of ~^qq% 
(see Fig. |^). Keeping at 100 GeV and increasing Y to 2, we find that for X = —2 
(corresponding to a Type I model with tan/5 = 0.5), the branching ratio is negative for 
all values of /ife G [2.4, 9.6] GeV. 

Going to heavier Higgs masses, e.g. mn = 500 GeV, we find for Y = 1 the same 
qualitative features as in Fig. ^, but shifted to larger ranges of |X|. We should warn, 
however, that for large enough |X|, when the Higgs contribution dominates over the 
SM one, the stability of the branching ratio becomes worse. The NLO corrections have 
a more dramatic dependence on the particular matching scale chosen, than in the case 
rUfj = 100 GeV. Indeed, for fi^^ = tt^h = 500 GeV, the relevant correction for the branching 
ratio 2Re{ADxY) is 70%, but exceeds -100% for = M^. 
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Figure 5: Branching Ratio as a function oiniu in a 2HDM of Type II, for three different choices of the 
low-scale: /ib = 2.4, 4.8, 9.6 GeV and /ivv = (upper frame), and for for three different choices of 
the matching scale: — M^, rrit, Trijj and /if, = 4.8 GeV. The needed input parameters are fixed at 
their central values listed in Table 1^ in Appendix P. 
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Figure 6: Contour plot in (tan P, TUh) obtained by using the NLO expression for the branching ratio 
BR(i? — >■ Xg'j) and possible experimental upper bounds (see text). The allowed region is below the 
corresponding curves. 

4.3. Type I and Type II Models 

Theoretical predictions for the branching ratio in Type II models stand, in general, on 
a rather solid ground. In Fig. |^, we show these predictions in a Type II model with 
tan/5 = 2, for 100 < ttlh < 600 GeV. This value of tan/? is particularly interesting since 
already for tan/3 > 2, the branching ratio becomes insensitive to the actual value of this 
variable: the contribution Dyy is, in fact, suppressed by the coupling |Yp < 1/4, whereas 
the contribution Dxy is multiplied by the coupling XY* = 1, for any value of tan/5. For 
tan/3 < 2, the branching ratio grows very rapidly when tan/? decreases and can be made 
compatible with existing measurements only for large values of 771^. 

The upper frame of Fig. ^ shows the low-scale dependence of BR{B Xg'j) for 
matching scale fiw = "^h^ 

for rriH > 100 GeV. It is less than 10% for any value of itih 
above the LEP lower bound of 45 GeV. Such a small scale uncertainty is a generic feature 
of Type II models and remains true for values of tan/? as small as 0.5. The lower frame 
in Fig. ^ shows the (very weak) matching scale dependence as obtained by varying /Xj^ 
in the interval [M^/, Max(mt, m^)]. We should point out that the lowest estimate of 
BK{B — i> Xs'~f) comes from the largest value of in this interval. We disagree on this 
point with ref. [|l7l where the lowest estimate is obtained with the smallest value of /i 



The result of []T^ is presumably due to the neglect of the term proportional to ln(yU^/M, 



w- 
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Figure 7: Branching Ratio as a function of nijj in a 2HDM of Type I, for three different choices of the 
low scale: jj^i, = 2.4, 4.8, 9.6 GeV and for /i^y = TTljj. The needed input parameters are fixed at their 
central values listed in Table |^ in Appendix ^. 

for the coefficients Cl''^^{fiw) and Cl''^^{fiw) in eq. (|17|). 

In Type II models, the theoretical estimate of BR(i? Xs'y) can be well above the 
range (1 — 4.2) x 10"'^ indicated by the CLEO Collaboration as the band of experimentally 
allowed values. It is therefore interesting to establish with some accuracy which values of 
{tan/3 , m^} are excluded by possible measurements of the branching ratio BR{B — > Xg'j) . 
Our results are given in Fig. ^ where we show the contour of the region excluded by the 
upper bound obtained by CLEO, 4.2 x 10~^, as well as for other two hypothetical values, 
3.4 X 10~^ and 5.0 x 10~^; the latter one is not far from the upper bound obtained 
by the ALEPH Collaboration. These contours are obtained by finding the minimum of 
BR(i? Xg^) , when varying simultaneously the input parameters within their errors 
(see Table |I] in Appendix |^) and the two scales /i^ and fiw in the ranges 2.4 < fib ^ 
9.6 GeV and M^i,. < fi^^ < Ma.x{mt,mH). For tan/? = 0.5,1,5, we exclude respectively 
'tTT'H ^ 375, 289, 255 GeV. Notice that the flatness of the curves shown in Fig. ^towards the 
higher end of m^, causes a high sensitivity of these bounds on all details of the calculation. 
Different treatments of the infrared singularities when the photon gets soft (e.g. with a cut 
in the photon energy |^) and the possible expansion of 1/Tsl in powers of could alter 



the branching ratio at the 1% level, i.e. well within the estimated theoretical uncertainty. 
These details, however, could still produce shifts of several tens of GeV, in either direction, 
in the lower bounds quoted above. 
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Branching ratios in Type I models can be reliably predicted for rriu > 100 GeV and 
tan/5 > 1. Theoretical results for this range of parameters fall within the CLEO band 
(1 - 4.2) X 10"^ as it can be seen in Fig. ^ for tan/5 = 1 and 100 < m„ < 600 GeV. 
Larger values of tan/3 decrease the Higgs contribution to the branching ratio, giving 
therefore values closer to the SM prediction. Lower values of rrijf (and/or tan/5 < 1) can 
produce results outside the CLEO range. For these parameters, however, the theoretical 
predictions are unstable under scale variation and, at times, ill-defined. As shown in 
Fig. 1^, a scale dependence of is obtained for itlh = 45 GeV and tan/5 = 1. For 
tan/5 < 1, negative values of the branching ratio are found already for = 100 GeV. 



4.4. Branching Ratio: Complex Couplings and Rate Asymme- 
tries 

In this section we study several aspects of the decay B — >• Xg'j in the presence of complex 
X and Y couplings. Particularly motivating is the observation that the measurement of 
BR(i? —>■ Xs'j) yields strong constraints on Im(XY*) pT] , p2| . In turn, these constraints 



limit the possibility of having large indirect CP asymmetries in neutral B decays such as 
B —y ipKg ||3^. This observation is based on a LO analysis whose reliability and stability 
under scale variation was not enquired. We plan to investigate this aspect and to check 



how the bound in [R^] may be modified by the inclusion of NLO corrections. 



The bound of ref. |^ on Im(XY*) is obtained as follows. Since the LO branching 
ratio is proportional to |C7''^^(/ife)p, the upper bound BR(i? X^'^) < 4.2 x 10~^ from the 
CLEO measurement implies an upper bound on \C^''^^{fif,)\'^, and therefore on Im(XY*) 
when: 

C?'^'^!/^.) = {c'rfM + |YpC°vt (/i.) + Re (X Y*)C?;^^(/i,)} + z Im (X Y*)C°;;^(/i,) 
= ^lm{XY*)C',fM, (45) 

i.e., when the real parts of the charged Higgs contributions cancel the SM coefficient 
Cj'gfj{fib). An inspection of the upper frame of Fig. ^ shows that, for rriu = 100 GeV 
and Y = 1, a vanishing branching ratio is induced by the the real coupling X = —2. 
The choice of complex couplings X and Y with Y = 1 and ReX = —2, therefore, fulfills 
to a good approximation the condition (p5[), for all values of ImX. The corresponding 
branching ratio, obtained with the central value of the input parameters, fiw = ^h, and 
Aife = f^b, and shown in Fig. |^ as a function of ImX, equals 4.2 x 10~^ at |ImX| ~ 2. For 
the chosen Higgs mass, mn = 100 GeV, therefore, the upper bound on |Im (X Y*)| is ~ 2. 
Notice that this procedure yields only a first estimate for the actual LO bound, since the 
errors of the input parameters have not been considered. 

We observe that a variation of the low-scale /if, in the usual range leads to large 
uncertainties for the branching ratio, throughout the whole range of Im X. In particular, 
for |ImX| ~ 2, this uncertainty amounts to ±25%. When including NLO corrections. 
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Figure 8: Branching Ratio as a function of IinX, with fixed ReX — — 2 and Y — 1, for — TTT'h 
and different low-scales: /ife = 4.8 GeV (solid lines), = 2.4 GeV (dashed lines) and yUj = 9.6 GeV 
(dash-dotted lines). The needed input parameters are fixed at their central values listed in Table ^ in 
Appendix^ Superimposed is the range of values allowed by the CLEO measurement. The upper (lower) 
frame shows the LO (NLO) result. 
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the situation does not improve, as the lower frame of Fig. ^ shows. It is interesting to 
see that the value of BR(i? — >■ Xg^) for iib = '"^b, at |ImX| ~ 2 drops from 4.2 x 10~^ 
to roughly 1 x 10~^. The intersection of the NLO curve for fif, = rrih with the CLEO 
upper bound is at ImX ~ 4. This procedure is essentially the construction of the NLO 
bound for |Im (XY*)| in the sense specified before. (In the actual construction, one should 
have cancelled also small real parts of C^'*^^ with NLO as and small terms coming from 
absorptive parts of virtual corrections.) 

The inclusion of NLO corrections verifies explicitly the instability of the LO upper 
bound on Im(XY*). Given the fact that the NLO predictions for BR(S Xs'y) are 
plagued by even larger scale uncertainties, it is hard to believe that the NLO candidate 
qualifies reliable bound. 

Not all complex X and Y couplings yield NLO predictions as problematic as those 
shown in Fig. A typical case in which the perturbative expansion of BR(i? — Xg'y) 
can be safely truncated at the NLO level is identified by: Y = 1/2, X = 2exp(20), and 
uIh = 100 GeV. Indeed, for these parameters, the real and imaginary part of \Y\'^Dyy and 

(xy*);d 

XY are dominated by Dgj^^. Therefore, the low— scale dependence of the branching 
ratio is not much larger than the very mild one obtained for the SM estimate, as shown in 
Fig. |. This case is particularly interesting since it gives rise to a theoretical prediction of 
BR(i? —>■ Xsj) consistent with the CLEO measurement, even for a relatively small value 
of rriH- Such a light charged Higgs can contribute to the decays of the t-quark, through 
the mode t —>■ H^b. 

We now investigate CP asymmetries induced by complex couplings. It is well known 
that in the SM a non-vanishing direct CP rate asymmetry 



is due to non-trivial relative weak- and relative strong phases in the decay amplitude for 
b —>■ s'-f (as well as the one for b — > s'-fg). We write the amplitude of 6 S7 and of the 
CP conjugated process as: 



where the dependence on the CKM matrix is manifest. Working in the LO approxima- 
tion, the quantities Au, and At are real, hence \A(b — >■ 57) p = \A{b 57) p, and 
consequently a^p = 0. Relative strong phases among Au, A^ and At only appear at the 
NLO level, due to absorptive terms in the loop diagrams. Thus, a non-vanishing CP 
asymmetry results if VubV*^, VcbV*^, and VtbVt* have also relative phases. In the approx- 
imation VubV*g = 0, which we used so far, the CP asymmetry vanishes, since unitarity 
implies that VcbV*^ = —VtbVt*. Undoing this approximation, the SM CP rate asymmetry 
turns out to be below 1% [@. 



BR{B Xs-f) - BR{B X,7) 
BR(B X,7) + BR{B -> Xsl) 



(46) 



A{b S7 ) 
A{b 57) 



(KbKTJ Au + (V^bV*) A, + (VtbVtl) At ; 

{VubV:,rAu + {v,bV:,rA, + (v,,v;jm, , 



(47) 
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Figure 9: Branching Ratio as a function of (f), where (j) parametrizes X: X = 2 exp(i0), for Y = 1/2, 
rrijj = 100 GeV, and fi^r = Tn^j. Sohd, dashed and dash-dotted hnes correspond respectively to 
/ifc = 4.8, 2.4 and 9.6 GeV. The needed input parameters are fixed at their central values listed in 
Table |l| in Appendix ^. Superimposed is the range of values allowed by the CLEO measurement. The 
upper (lower) frame shows the LO (NLO) result. 
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Figure 10: CP rate asymmetry ttcp as a function of 0, where parametrizes X: X = 2exp('j0), 
for Y = 1/2, rriu = 100 GeV, and = rriu. Solid, dashed and dash-dotted Unes correspond 
respectively to fib = 4.8, 2.4 and 9.6 GeV. The needed input parameters are fixed at their central 
values listed in Table |^ in Appendix In the upper (lower) frame the LO (NLO) expression for the 
denominator in eq. (EgI) is used. See text for more details. 
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Since the CP rate asymmetry is so tiny in the SM, it would be exciting if the imaginary 
parts of the X and Y couphngs would induce -together with the absorptive parts of 
the NLO loop-functions- measurable CP rate asymmetries. In order to investigate this 
question, we switch off the SM asymmetry by working in the limit VubV*^ = 0. As an 
illustrative example we calculate the CP asymmetries for the same values of parameters 
as in Fig. ^, i.e., for X = 2exp(20), Y = 1/2 and for = 100 GeV, where the branching 
ratio appears to be well-behaved. As mentioned above, the NLO prediction for the 
branching ratio in the numerator of eq. ( ^61) is required to obtain the first non- vanishing 
term for the rate asymmetry. For the denominator, one can either use the LO or the 
NLO estimate for the two branching ratios. The difference amounts to higher order terms 
which are not systematically calculated. The respective results are shown in the upper 
and lower frame in Fig. |10| as a function of 0. For this specific choice of X and Y the 
asymmetries are rather modest, at the 1% level. The scale dependence shown by the lower 
frame is certainly smaller than that in the upper frame. However, since both procedure 
are in principle viable, one cannot conclude that the result for the CP asymmetries is as 
reliable as indicated in the lower frame. Although we have not systematically scanned 
the parameter space {X, Y,m^}, this is our generic result: choices of couplings X and Y 
which render the branching ratio stable, i.e. not plagued by a large dependence on the 
scale fj,b, induce small values of acp. 



5. Conclusions 



We have presented NLO predictions for the decay B Xgj in generic 2HDMs with 
flavour-conserving tree-level neutral couplings. We include in this definition Multi-Higgs 
Doublet models where all charged Higgs bosons, except one, are heavy and completely 
decoupled at the electroweak scale. Their existence leaves an imprint only in the Yukawa 
potential where the couplings X and Y multiplying respectively the down- and up-term 
are not necessarily correlated, and in general complex. This generalization allows a si- 
multaneous study of -B — s> for different types of 2HDMs by continuously varying the 
couplings X and Y. Results for the well-known Type I and Type II models are then 
obtained for specific choices of real couplings: only one combination of X and Y, usually 
denoted as tan/3 , occurs. Since supersymmetric models have an enlarged Higgs sector 
with two Higgs doublets of Type II, the results presented here are also a first step towards 
a complete evaluation of the rate of B Xsj in these models. They constitute already a 
good approximation for those supersymmetric scenarios, such as gauge-mediated super- 
symmetric models with large tan/? , where the Higgs contribution dominates by far over 



the genuinely supersymmetric contributions 



Our calculation is carried out using the effective Hamiltonian formalism with on-shell 
operators. The NLO matching condition for the decay B —>■ Xg'j was already completely 
calculated in and partially in |3^. For the Wilson coefficients Cj^{fiy^,) and C^^{fivy), 



we find agreement with the existing results, but we disagree with ref. |27,17] for the 



28 



matching scale dependence of the coefficients Cf*^(/ivv) and C^^^fi^)- We generahze the 
solution of the RGE needed to obtain the coefficient C^^^Hb) at the low-scale Hb for values 
of 

fiw different from M^. We also correct the dependence on the Wilson coefficients of 
one of the non-perturbative contribution ioT[B ^ ^sl)-, which is erroneously reported 

m 



m 



We have given predictions for BR(i? Xg'y) in 2HDMs as functions of the param- 
eters {X^Y^uih}- We found that the theoretical uncertainties of these NLO calculations 
are in general larger than those obtained in the SM. The quality of our predictions, there- 
fore, depends strongly on the values of the parameters considered. As a ffist step, we 
summarize the theoretical features of our results, ignoring a comparison with the existing 
experimental data. We distinguish several cases: i) When these parameters are such that 
the SM contribution is much larger than the Higgs contribution to BR(-B Xg'-)) , the 
reliability of our predictions does not differ much from that in the SM. This, however, 
does not prevent the BR(i? — > Xs7) from being rather sensitive to m^. ii) For values of 
parameters which bring the Higgs contributions to the same level as the SM one, with 
constructive interference between the two, we find in general larger scale dependence than 
in the SM. iii) When the Higgs contributions are still of the same order as the SM con- 
tribution, but interfere destructively with it, we find very large scale dependences and, 
in specific cases, we obtain negative values of BR(i? — > Xg'-y) . iv) When X and Y make 
the Higgs contributions largely dominate over the SM one, we find a scale dependence in 
BR(5 — >• Xg'j) of ~ 40% for mjj of C(Mj^). For much larger values of m„, the quality 
of the theoretical prediction depends strongly on the matching scale considered. 

Type II models are typical of case i) and barely become of case ii) for itIh at the 
LEP lower bound of 45 GeV |^ and for tan/? ~ 0.5. The theoretical uncertainty for such 



models is in general below 10%. Similar uncertainties are obtained for Type I models 
when rriH and/or tan/? are large enough. Rather unstable results are, however, obtained 
for ttIh at the lower end of allowed values and tan/3 ^ 

When these problematic regions are avoided, then all predictions obtained for Type I 
models are consistent with the CLEO measurement of BR(S — > Xg'j) . The same is not 
true for Type II models where the theoretical predictions are always above the SM result. 
Measurements become then highly constraining. Taking into account only the CLEO 
result, without combining it with the still preliminary one from ALEPH, we exclude 
TTT'H < 375 GeV for tan/5 = 0.5 and itih < 255 GeV for tan/3 = 5. These bounds are very 
sensitive to details in the definition of the branching ratio (e.g. whether it requires or not 
an expansion in ag of l/Fg^), which give rise to uncertainties usually not included in the 
theoretical error for this observable. For the generic 2HDMs which we consider, we find 
wide regions of parameter space where the theoretical predictions for BR(i? —>■ Xgj) are 
reliable and within the band of values (1.0—4.2) x 10~^ allowed by the CLEO measurement. 
In particular, we find that within these models, charged Higgs bosons can still be light 
enough to be produced through a decay of the t-quark. 

Since X and Y are in general complex, new CP violating effects are induced. In partic- 
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ular, the combination Im (X Y*) is important as it may affect the indirect CP asymmetries 
in i? — > ipKg aheady at the LO. It is known that the measurement of BR(i? — Xg'^) con- 
strains Im(XY*) 1^ and an upper bound for this quantity has been obtained in 



using the LO calculation. We find that the LO estimate of the branching ratio at the 



values of X and Y which determine the upper bound extracted in is unstable under 
scale variation. The addition of NLO corrections tend to shift the LO bound to a higher 
value. Nevertheless, the scale dependence of BR(i? — > Xg^) for the new combination of 
X and Y needed, is still too large to conclude that the bound obtained at the NLO level 
is stable against higher order corrections. The large instabilities encountered even at this 
order are due to the fact that the construction of this upper bound requires the almost 
complete cancellation of the SM contribution. 

At the NLO, complex values of X and Y induce also direct CP rate asymmetries, 
in — > , which is essentially vanishing in the SM. If sizable, a measurement of this 
observable, could provide a handle to detect some of these models. Unfortunately, we find 
that combinations of couplings where BR(i? — > Xs7) is reliably predicted, lead to rate 
asymmetries only at the 1% level. 

We find that within 2HDMs, the truncation of the perturbative series at the NLO 
level is often inappropriate. This is in sharp contrast with the SM case, where the large LO 
theoretical uncertainties are drastically reduced in the NLO calculation, and the overall 
size of the NLO corrections is, in comparison, rather modest. It is somehow disturbing 
to find the problematic features described in this paper in models which structurally do 
not differ very much from the SM and it is conceivable that other extensions may suffer 
from the same diseases. One may wonder about the deeper reasons for the cancellations 
observed in the SM. Their understanding may help in clarifying also the reasons for the 
features encountered in 2HDMs. 
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Appendix A. Wilson Coefficients at /i^F 



We list here the functions introduced in the text, which define the Wilson coefficients at 
the matching scale (see eqs. (p!^)-([T8|)). 

SM case 

The LO functions are (x = m^/M^): 

-8x^ + 3x^ + 12X-7+ (18x2 - 12x) Inx 



X 

24 

X 



{x - 1)4 

-x^ + Gx^ — 3x — 2 — 6x In X 
(x- 1)4 



(Al) 



In the MS scheme, we have at the NLO: 



En 



8,SM 



M' 



7,SM 



x(x2 + llx-18) x2(4x2 - 16x + 15) , 2, 2 

-, r:^ 1 -, n In X In X 

12(x-l)^ 6(x-l)4 3 3 



-16x4 - 122x3 + 80x2 _ g 



X _ A 1 \ 6x4 + 4g^3 _ 28a;2 ^ ^ 



Li2 1 

X 



9(x-l)4 "V xj 3(x-l)5 

-102x5 - 588x4 _ 2262x3 + 3244x2 - 1364x + 208 , 

H ; In X 

81(x-l)5 

1646x4 + 12205x3 - 10740x2 + 2509x - 436 
^ 486(x- 1)4 

-4x4 + 40x3 + 41x2 + X / I 
Li2 1- 



In X 



6(x- 1)^ 



X 



-17x3 - 31x2 
2(x - 1)5 



In X 



-210x5 + 1086x4 + 4893x3 + 2857x2 - 1994x + 280 , 

H ; — mx 

216(x - 1)5 

737x4 _ 14102x3 - 28209x2 + 610x - 508 

^ 1296(x - 1)4 

82x5 + 301x4 + 703x3-2197x2 + 1319x-208 -(162x4 + 1242x3-756x2) Inx 

81(x- 1)5 

77x5 _ 475^,4 _ 11112.3 ^ go7a;2 ^ 1042^ _ 140 +(918x3 + 1674x2) Inx 

108(x- 1)5 

X [47x3 - 63x2 _^ 93, _^ 7 _ (ig2;3 _)_ 302;2 _ 24a;) In x 

(x - 1)5 
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T8,SM — 2x 



-x^ -9x'^ + 9x + l + (6x^ + 6x) \nx 
{x - 1)5 



(A2) 



2HDM case, coupling |yp 

The functions relative to charged Higgs contributions, with coupling |y|^, are at the LO 
(y = m2/m|): 



(A3) 



at the NLO: 



E 



H 



W; 



7,YY 



M' 



36 



9^ 



7y^ - 36y2 + 45y - 16 + (18y - 12) Iny' 



Li2 (1 ) H ^ — In y 



(y - 



y, 



(y - 1)^ 



21/ - 192y=^ - 174y2 + 251y - 50 



9{y - 1)5 
-1202y=^ + 7569y2 - 5436y + 797 



— Eff 
9 



7,r^ - 27^ 



8,YY — 



36^ 



108(?/ - 1)4 

42y4 + 318y=^ + 1353^/2 + 817y - 226 , 
+ 36fe3I)^ 

-4451y=^ + 7650y2 - 18153y + 1130 
+ 216(y - 1)4 

-142/4 + i49|/3 - I53y2 - iSy + 31 - {18y^ + 138^^ - 8Ay) Iny 
-7y^ + 25y^ - 279y'^ + 223?/ + 38 + (102^^ + I86y) Iny 



— 7 

6 



T, 



7,YY 



8,YY 



2 T7,sm{x 



y) 
y) 



(A4) 
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2HDM case, coupling (XY*) 

Similarly, the functions relative to the charged Higgs contributions proportional to (X Y*) 
are at the LO: 



12 
1 



,XY 



-5y'^ + - 3 + (6y - 4) In y 
-y"^ + % - 3 - 21n|/ 



(A5) 



and at the NLO: 



- V 
3^ 



Li2 (1 1 + 7T:^ln y 



3(1/ -1)= 



4y3-2V + 2i/ + 6 ^ 

H ^ n In y + 

3(y - 1)4 ^ 



-2|/^ + 13|/ - 7 ' 



M' 



8,xy 



7,XY 



14y3 - 12y2 + I87y + 3 3(2V - 44i/ + 143) 



7,xY = 77 y 



3^ 



4(y - 1)4 ^ 8(y - 1)3 

-8y^ + 55y2 - 68y + 21 - (Gy^ + 28y - 16) Iny 

-7y^ + 23y^ - 97y + 81 + {34y + 38) Iny 

13y^ - 20y + 7 - (6y^ + 4y - 4) In y^ 

{y - 1)^ 



8,XY 



= 2t/ 



-y^-4y + 5 + (4y + 2)ln|/ 



(A6) 



Appendix B. Anomalous Dimension Matrix 



For completeness, we give the anomalous dimension matrix which govern the evolution of 
the Wilson Coefficients from /i^ to /ib. It can be expanded perturbatively as: 



eff/ \ _ f^^sip} 0,eS , Q'KA'-) ^ l,eff 



7r(/x) 



47r 



+ 



(47r) 



7r + 0{af). 



2 o 



(Bl) 
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The matrix 7^^'^''^ is given by: 



-4 


8 
3 





2 
9 








208 
243 


173 " 
162 


12 








4 
3 








416 
81 


70 
27 











52 
3 





2 


176 

81 


14 
27 








40 
9 


100 
9 


4 
9 


5 
6 


152 
243 


587 
162 











256 
3 





20 


6272 
81 


6596 
27 








256 
9 


56 
9 


40 
9 


2 

3 


4624 
243 


4772 
81 




















32 
3 























32 
Q 


28 



and in the MS scheme with fully anticommuting 75, 7^^^ is 



(B2) 



355 


502 


1412 


1369 


134 


35 


818 


3779 
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27 


243 


243 


243 


162 


243 
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508 
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3 
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81 


81 


81 


27 


81 
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4468 


31469 
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3373 


22348 


10178 


81 


81 


81 
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243 
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8158 
243 


59399 
243 


269 
486 


12899 
648 


17584 
243 


172471 
648 








251680 
81 


128648 
81 


23836 
81 


6106 
27 


1183696 
729 


2901296 
243 








58640 


26348 


14324 


2551 


2480344 


3296257 


243 


243 


243 


162 


2187 


729 




















4688 
27 























2192 

SI 


4063 
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(B3) 



Appendix C. "Running" numbers 



The vectors {oi}, {hi}, and {a'j}, {h'^, needed for the evaluation of the low-scale Wilson 
coefficients CY'^^{fib) and Cg'^^dxb) are: 



{h^} 



those needed for C7' (nb): 



14 16 6 12 1 
23' 23' 23'~23' 0-4086,-0.4230,-0.8994, 0.1456 j 

626126 _56281 _3 _^ _o.6494, -0.0380, -0.0186, -0.0057 \ , (CI) 
272277 51730' 7 14' ' ' ' r ' v ; 

^, 0.4086, -0.4230, -0.8994, 0.1456 
3 13063 

— — , -0.9135, 0.0873, -0.0571, 0.0209 \ ; (C2) 
363036 

l.cff/ 
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, , r 4661194 8516 1 

{a} = , , 0, 0, -1.9043, -0.1008, 0.1216, 0.0183 \ 

^ ^ I 816831 ' 2217' J 

{/.} = {-17.3023, 8.5027, 4.5508, 0.7519, 2.0040, 0.7476, -0.5385, 0.0914} 

{ki} = { 9.9372, -7.4878, 1.2688, -0.2925, -2.2923, -0.1461, 0.1239, 0.0812} 

{/.} = { 0.5784, -0.3921, -0.1429, 0.0476, -0.1275, 0.0317, 0.0078, -0.0031} . (C3) 



Appendix D. Virtual correction functions 



The renormalization scale independent parts of the virtual corrections, encoded in the 
functions in eq. (pOD, read: 



r2 



243 



{-833 + 1447rV/^ 



+ 
+ 
+ 



1728 - ISO-K^ - 1296 C(3) + (1296 - 32A7t^)L + 108^^ + 36L' 



648 + 727r2 + (432 - 216n^)L + 36L' 



-54 - 847r^ + 1092L - 756L^ 



\ -5 + 

81 ^ 



} 



45 - 37c' + 9L 
+ 0{z') 



9L' 



z + 



+ [28 - 12L] } 



= — 



32 
9" 



-vr 
9 



4 



— ^(-33 + 27r2 
27^ 



Gin) 



(Dl) 



where z is defined as 2; = ml /ml and the symbol L denotes L = ln{z). Notice that 
^3 ,f4 ,f5, and tq are not used in the approximation C°''^^(/ife) =0 (i = 3, 4, 5, 6) used for 
the matrix elements. 



Appendix E. Bremsstrahlung fij terms 

The expressions for the bremsstrahlung functions fij we use in the present paper are 
obtained after integrating one variable in the expressions given in Appendix B in ref. 
TB|. The exphcit expressions read (converted to the operator basis |T2p: 



/ll — tjt: /22 
36 



35 



/l2 
/l7 
/l8 

/22 

/27 
/28 

/78 

/ss 



-^/22 

— 7 /27 


^7 /28 
O 

162 /-lA 



27 Jo 



dt{l - ztf 



mi 
t 2 



9 



I dtil - zt) [Git) + - 



1 


/27 






^25 


-i) 


-I 

9 ' 






1 


/16 




27 




3 



+ 41n 



(El) 



where z is z = ml /ml and /Xfe is the renormahzation scale. The function G{t) appearing 
in eq. (ETI) reads 



G{t) 



—2 arctan^ J ^ 



-7rV2 + 21n2 (V^±^) _2mln( 



for t < 4 
for t > 4. 



(E2) 



Appendix F. Input Parameters 

We hst in table |I] the values of input parameters used in our calculation. The masses 





mt 
(GeV) 


mc/mb 


rrib - mc 
(GeV) 






BR,L 


0.119 


175 


0.29 


3.39 


130.3 


0.95 


0.1049 


±0.004 


±5.0 


±0.02 


±0.04 


±2.3 


±0.03 


±0.0046 



Table 1: Central value of our input parameters (first line) and their uncertainties (second line). 



mj, m^, and are understood to be the pole masses of the top, bottom and charm 



quark. The value of the mt is obtained by combining the results given in ||50| , pl| ; those 
of the two combinations of and m^ are taken from [Q. For ctem and |l^*VJf,/V^feP 



we refer to |^ and respectively. We take as(M^) = 0.119 ± 0.004, as an average 
between a pessimistic and optimistic estimate of the error [54|, as suggested in |55|. The 
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value of the semileptonic branching ratio BRc;^ has been recently obtained by the CLEO 
Collaboration |^. The other constants used in the calculation are = 80.33 GeV, 
Ai = -0.5 GeV^, and A2 = 0.12 GeV^. 
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